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Abstract 

We describe an efficient position space technique to calculate lattice Feyn¬ 
man integrals in infinite volume. The method applies to diagrams with mass¬ 
less propagators. For illustration a set of two-loop integrals is worked out 
explicitly. An important ingredient is an observation of Voliwinkel that the 
free lattice propagator can be evaluated recursively and is expressible as a 
linear function of its values near the origin. 
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1. Introduction 


Perturbation theory plays an important conceptual and practical role in 
lattice field theory. The continuum limit can be studied in this framework and 
numerical simulations must often be supplemented by perturbative calcula¬ 
tions to obtain the desired physical result. The evaluation of lattice Feynman 
diagrams is not easy in general, because the Feynman integrands tend to be 
rather complicated functions of the loop and external momenta. Standard 
tools, such as the Feynman parameter representation and partial integration 
methods, are usually not very helpful on the lattice. 

In a recent computation of the relation between the bare lattice coupling 
and the MS coupling in SU(iV) gauge theories [1] we have made extensive 
use of some new position space techniques. They are expected to be more 
generally applicable and so it is our aim in the present paper to describe these 
methods in some detail. Basically we discuss a number of examples and leave 
it to the reader to adapt the methods to the case of interest. In other words, 
we do not attempt to set up a general strategy for the computation of lattice 
Feynman diagrams. We feel that more experience must be accumulated before 
this is possible. 

Throughout this paper we assume that the theory considered lives on a 
four-dimensional hypercubic lattice A with spacing a = 1. We only discuss 
diagrams with free propagators equal to 1/p 2 , where p = (po,.. • ,P 3 ) denotes 
the line momentum and 


3 

Pn = 2 sin {\p^) , P 2 = YjPI- 

/i=0 


( 1 . 1 ) 


This covers the case of (pure) lattice gauge theories with standard Wilson 
action. Fermion propagators are different and some additional work will be 
required before our methods can be applied to lattice QCD. 

The basic ideas to be discussed in this paper are best explained by con¬ 
sidering the simple two-loop integral 


p d 4 k d 4 q 1 
(2vr) 4 (2 tt) 4 k 2 q 2 r 2 ’ 


( 1 . 2 ) 


Aq can be calculated directly using a Monte Carlo program, for example, but 
it would be difficult to obtain the integral very accurately in this way. Now if 
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we introduce the position space propagator 


. f* d 4 p e ipx 
J-* (2tt) 4 p 2 

the integral may be rewritten in the form 

A 0 = ^ G{xf. 

xeA 


(1.3) 


(1.4) 


Since G(x) falls off like 1/x 2 at large x, the sum is absolutely convergent and 
a first approximation to Aq may be obtained by summing over all points x 
with say \x\ < 20. Taking the detailed asymptotic form of G{x) into account, 
one can then systematically improve on the accuracy until the desired level of 
precision is reached. 

For this procedure to work out one must be able to calculate the prop¬ 
agator G{x) accurately at all lattice points x that are not too far away from 
the origin. An elegant solution to this problem has been given by Vohwinkel 
[2]. His method allows us to obtain the propagator to essentially any desired 
precision with a minimal amount of computer time. 

In sect. 2 we discuss the free propagator in detail and explain Vohwinkel’s 
method. We then address the question of how to calculate lattice sums accu¬ 
rately and consider various examples to illustrate the procedure (sect. 3). The 
position space method is also useful to calculate Feynman diagrams with non¬ 
zero external momenta, particularly in the continuum limit. This is demon¬ 
strated in a number of cases in sect. 4. Some concluding remarks are collected 
in sect. 5. Much of the technical material has been deferred to a series of 
appendices to make the paper more readable. 


2. The free propagator 


We begin by discussing the behaviour of G{x) at large x and then derive 
a recursion relation which allows us to express the propagator through its 
values close to the origin. After some further steps a strategy for the numerical 
calculation of G{x) emerges which proves to be very efficient and simple to 
program. 
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2.1 Asymptotic expansion at large x 

In terms of the forward and backward lattice derivatives, 


V^fix) = f(x + p) - f(x), Vlf(x) = f(x)-f(x-p) (2.1) 

(where fi denotes the unit vector in direction p), the lattice laplacian is given 

by 

3 

A = ^V*V,. (2.2) 

/i=0 

G{x) is a Green function for A, i.e. it satisfies 

-AG(x) = { 1 lfa; = 0 ’ (2.3) 

L 0 otherwise. 

In the continuum limit (which here amounts to taking x to infinity) we thus 
expect that G{x) converges to (47r 2 a; 2 ) -1 , the Green function of the continuum 
laplacian on 1R 4 . 

A rigorous derivation of the large x behaviour of G{x) confirms this (ap¬ 
pendix A). One may also work out the subleading terms and obtains the 
asymptotic expansion 


G(x) 


4ir 2 c 


1 x 4 

s 2 ^ (x 2 ) 3 


„ 4\2 


- 4- 


,,232 


16- 


,234 


- 48- 


,235 


40- 


,236 


(2.4) 


which includes all terms up to order \x\ 6 . For notational convenience the 
shorthand 


x 


n 




l-i=0 


(2.5) 


has been used in eq.(2.4). 
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2.2 Recursion relation 

The key to Vohwinkel’s method is the observation that 

(V$ + V ll )G(x) = x ll H(x), (2.6) 

where 

H( x) = J^0 } -^Hf) (2.7) 

is independent of p. Eq.(2.6) thus represents a set of 4 relations for G and H 
of which one may be used to eliminate H. For example, by summing over p 
and using eq.(2.3) one obtains 

3 3 

H(x) = (2/p)^2 [G(x) ~ G(x - p)] , p=^x ll , (2.8) 

fi =0 fi =0 

and hence 

3 

G(x + p) = G(x — fi ) + (2 x^/p) [G(x) — G(x — />)]. (2-9) 

v=0 

Of course this relation is only valid at the lattice points x where p does not 
vanish. 

The propagator G{x) is independent of the sign and the order of the 
coordinates (xq, ..., x^). We can hence restrict attention to the points x with 
x o > x \ > X 2 > S 3 > 0. In this sector eq.(2.9) is a recursion relation which 
allows us to express G{x) as a linear combination of 

0 ( 0 , 0 , 0 , 0 ), 0 ( 1 , 0 , 0 , 0 ),..., 0 ( 1 , 1 , 1 , 1 ), ( 2 . 10 ) 

i.e. by the values of the propagator at the corners of the unit hypercube. 

The recursion alone does not imply any restrictions on these 5 initial 
values. When further properties of the propagator are taken into account one 
may however show that 


G(0,0,0,0)- G(l, 0,0,0) = 1/8, (2.11) 

G(0, 0, 0,0) - 3G(1,1, 0,0) - 2G(1,1,1,0)= 1/tt 2 , (2.12) 

G(0,0,0,0)- 6G(1,1,0,0) - 8G(1,1,1,0) - 3G(1,1,1,1) = 0. (2.13) 
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The first of these relations is equivalent to eq.(2.3) at x = 0. An amusing way 
to deduce the other two is given in appendix B. Thus, from the 5 initial values 
3 can be eliminated and the propagator is obtained in the form 

G(x) = ri(a;)G'(0, 0, 0,0) + r 2 (x)G( 1,1,0,0) + r 3 (x)/ir 2 + 7 - 4 ( 2 :), (2.14) 

where 7 - 1 ( 2 :),..., 7 - 4 ( 2 :) are rational numbers recursively computable through 
eq.(2.9). 

2.3 Numerical computation of G{x) 

Using MAPLE or some other algebraic programming language, it is straight¬ 
forward to write a program which calculates the coefficients 7 - 1 ( 2 :),.. ., 7 - 4 ( 2 :) 
exactly. We are then left with the task to compute G(0, 0, 0,0) and G{ 1,1, 0,0) 
numerically. It is necessary to obtain these numbers very accurately, since 
the significance losses in eq.(2.14) can be huge. The reason for this is that 
7 - 1 ( 2 :),..., 7 - 4 ( 2 :) are exponentially increasing at large x while G{x) remains 
bounded. 

An interesting observation now is that we can turn the tables and use the 
rapid growth of the coefficients 7-1(2:),..., 7-4(2:) to calculate G(0, 0,0,0) and 
G{ 1,1,0,0) in a most efficient way. To this end we consider the lattice points 
x = (n, 0, 0, 0) and y = (n, 1, 0, 0), where n is an adjustable integer. Inspection 
shows that the associated coefficients rk{x) and rk(y) are roughly proportional 
to 10”. So if we express G(0, 0, 0,0) and G{ 1,1, 0,0) as linear combinations of 
G(x), G(y ), 1 /7r 2 and 1 [by solving eq.(2.14) and the corresponding equation 
for G(y)\, the coefficients multiplying G{x) and G(y) will be of order 10“”. 
Neglecting these terms one thus obtains G(0, 0,0,0) and G(l, 1,0,0) to this 
precision. 

Note that this method of calculating G(0, 0,0,0) and G(l, 1,0,0) is ex¬ 
ponentially convergent. Moreover the error can be estimated rigorously since 
the neglected terms are known and can be bounded easily. With a negligible 
amount of computer time one thus obtains 

G(0,0,0,0) = 0.154933390231060214084837208107 ..., (2.15) 

G( 1,1,0,0) = 0.012714703770934215428228758391..., (2.16) 

and it would cost little more time to get the next 30 digits of these numbers. 
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3. Evaluation of lattice sums 


We now discuss how to calculate lattice sums such as the one appear¬ 
ing in eq.(1.4). As a preparation we first consider the zeta functions Z(s,h) 
associated with the hypercubic lattice A. 

3.1 Zeta functions 

Let h(x) be a harmonic homogeneous polynomial in the real variables x = 
(Aq, ..., X 3 ) with even degree d. Examples of such polynomials are 


h 0 (x) = 1, (3.1) 

hi(x) = 2a ; 4 — (a; 2 ) 2 , (3-2) 

h 2 (a;) = 16a; 6 — 20a; 2 a ; 4 + 5(a; 2 ) 3 , (3-3) 

h 3 (x) = 560(a; 4 ) 2 - 560a; 2 a; 6 + 60(a; 2 ) 2 a; 4 - 9(a ; 2 ) 4 (3.4) 

[cf. eq.(2.5)]. For complex s with Res > 2 + d /2 we dehne a generalized zeta 
function through 

Z(s,h) = Yj h{x){x 2 )-% (3.5) 

xeA 

where the primed summation symbol implies that the point x = 0 should be 


omitted. Note that the sum (3.5) is absolutely convergent in the specified 
range of s. 

The zeta functions can be calculated as follows. We first introduce the 
heat kernel 

k(t,h) = h(x)e~ wtx (3-6) 

xeA 

and rewrite eq.(3.5) in the form 

7T S f°° 

z (s,h)=-— dtt 8 - 1 [k(t,h) - h(0)]. (3.7) 

1 Jo 

Using the Poisson summation formula 

^2 e~ iqx e~' Ktx2 = t ~ 2 e -(a+ 2 ^) 2 /(4 *t) 

x£A x£A 
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and the property that h is harmonic, we deduce that the heat kernel obeys 

k(t,h) = (-l) d/2 t- d - 2 k(l/t,h). (3.9) 

Now if we split the integral (3.7) in two parts and apply the identity (3.9), the 
representation 


Z(s, h) 


7r s f 2h(0) 

rR1^2) 


+ dtit 8 - 1 + (-l) d / 2 t d - s+1 ][k(t,h)-h(0)]^ (3.10) 


is obtained. In the integration range 1 < / < oo the series (3.6) is rapidly con¬ 
vergent. Eq.(3.10) is hence suitable for numerical evaluation. In the cases of 
interest the zeta functions can in this way be calculated to (say) 16 significant 
digits with little effort. 


3.2 Computation of Aq 

Since the propagator G{x) is known accurately, we can obtain an estimate of 
the sum (1.4) by summing all terms with \x\ < 20. The result of this calcu¬ 
lation is accurate to a relative precision of about 10 —4 , i.e. the approximation 
is already quite good. One could of course do better by including more terms, 
but since the sum is only quadratically convergent this is not the most efficient 
way to proceed. 

The convergence can be accelerated by making use of the known asymp¬ 
totic form of the propagator. From eq.(2.4) we deduce that 


G(xf = {GOrflas + 0(M- 12 ), 

x —> OO 


(3.11) 


where 




(47T 


233 


,,233 


10(x 


235 


h 0 (x) 


24 


,,2 36 


7(x 


2)7 


hi(x) 


3 

A{x 2 f 


h 2 (x) 


33 

140(a; 2 ) 9 



The subtracted sum 


A' = [G(x f - 

xeA 


(3.12) 


(3.13) 
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is hence rapidly convergent and so may be computed accurately by summing 
over all points x with \x\ < 20. The subtraction term 

.4;' = ]T'{G(.tf}„ (3.14) 

xeA 


may be expressed as a sum of zeta functions, 

AS = ("4 7 r 2^3 ho) ”1” Yo^ 5, ^ ^i) 

24 S qq 'I 

+ y Z(7,M - 4^(8,M + —z(9,h 3 )|, (3.15) 

and is thus accurately calculable, too. Noting 

Ao = A' 0 A Aq + Gr(0) 3 , (3.16) 

the final result then is 

A 0 = 0.0040430548122(3). (3.17) 


A more precise determination of Aq along these lines would of course be pos¬ 
sible, but it seems unlikely that this will ever be required. 


3.3 More complicated integrals 

The techniques which have been discussed so far can easily be generalized to 
compute integrals of the form 


d 4 k d 4 q P(k,q,r) 

-x (2vr) 4 (2 tt) 4 


(3.18) 


where P(k , q, r) is a polynomial in the sines and cosines of the momenta k, q, r 
and n,m,l are positive integers (r = —k — q as before). To illustrate this we 
consider the integrals 


Q i 


j"* d 4 £; d 4 g 

-*■ (2vr) 4 (2 tt) 4 ^WqH*' 


(3.19) 
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<£> 

to 

II 

r d 4 k d 4 q A 

~2 ~2 
hj /i L lli l /I 

(3.20) 

U (2vr) 4 (2vr) 4 E 

k 2 q 2 r 2 

Ri = 

r d 4 k d 4 q A 

~2 

^li^/i 

(3.21) 

U (2vr) 4 (2vr) 4 E 

(k 2 ) 2 q 2 r 2 ’ 

ii 

r d 4 k d 4 q A 

K 

(3.22) 

U (27T ) 4 (2vr) 4 E 

{k 2 fq 2 r 2 


In position space the first of these is given by 

3 

Qi = E E [-' KV,G(x)f G(x) (3.23) 

x£A /i=0 

[cf. eq.(2.1)]. The terms in this sum fall off like |a;| —10 and a direct evaluation 
without subtraction yields 

Qi = 0.0423063684(1). (3.24) 

In the same way the result 

Q 2 = 0.054623978180(1) (3.25) 


is obtained. 

To compute the integral (3.21) we introduce the auxiliary function 

A(l) = /' (sy T" _ 9 / (#2)2 - (3 - 

The position space representation of R\ then reads 

3 

ifi = EE [-V*V At ^r(a:)] [~V$V»G(x)] G(x). (3.27) 

x£A ii=0 

Some properties of the function K(x) are discussed in appendix C. In par¬ 
ticular, it is shown there that K is recursively computable in terms of the 
propagator G. We may thus apply the techniques explained above and obtain 

R x = 0.006603075727(1), (3.28) 
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where a subtraction has been made to accelerate the convergence of the sum 
(3.27). 

In the case of the integral (3.22) another auxiliary function, 

L( x ) = f (e vpx - 1 - ipx + \{pxf) /(p 2 ) 3 , (3.29) 

must be introduced. The notation 

3 

p ll = sin(p ll ), px=^ j p ll x ll , (3.30) 

/i=0 

has here been used. This function is also calculable in much the same way as 
K (appendix C). In position space we have 

3 

T 3 = E E [VSV.VJV^*)] G(x)\ (3.31) 

x £A /j ,=0 

and the result 

T 3 = 0.00173459425(1) (3.32) 

is now obtained straightforwardly by applying the subtraction method. 


4. Feynman integrals with non-zero external momenta 


In this section we consider the Feynman diagrams shown in fig. 1. The 
external momentum is denoted by p and the numerators of the corresponding 
Feynman integrands are set equal to 1. We are interested in the continuum 
limit of these integrals. Since we are using lattice units and since there are 
no mass parameters, the continuum limit is equivalent to taking p to zero, 
i.e. our task is to work out the leading terms in an asymptotic expansion of 
the diagrams for p —► 0. 


10 



—O— Q ^3^- 

a b c 

Fig. 1. Self-energy diagrams studied in sect. 4. (b) and (c) are referred 
to as the “bigmac” and the “eye” diagram, respectively. 

4-1. Scalar one-loop integral 

The integral corresponding to the first diagram in fig. 1 is given by 

, r d4fc 1 

Mf) = Lw?W’ s = k + p ' (4 ' 11 

A(p) is a smooth function at all momenta p that do not vanish modulo '2ir. 
We would like to show that 


A{p) pZo (4^)2 [— ln (l j2 ) + 2 ] + p 2 + 0(p 2 ), (4.2) 


where Pi is a constant [in this section 0 (p n ) stands for a remainder R(p) such 
that linip^o R(p)/\p\ n ~ e = 0 for all < > 0]. 

We first compute A(p) in a conventional way using the Reisz power count¬ 
ing theorem [3,4]. The idea is to subtract the integral at p = 0 to make it 
convergent in the continuum limit. We thus define a subtracted integral, 



1 

k 2 s 2 


1 


( k 2 + m 2 ) 


2 


1 


and the subtraction term 



d 4 k 1 

(2vr) 4 (jfc2 + m 2)2’ 


where m > 0 is an auxiliary mass parameter. The full integral 


A(p) = A'(p, m) + A"(m) 


(4.3) 


(4.4) 


(4.5) 


is of course independent of m. We are thus free scale m to zero proportionally 
to p. I 11 this limit the Reisz power counting theorem applies and one concludes 
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that 


Urn A\\p,\m) = j_^ AA {^ ~ 

= jhf n ™ 2 / p 2 ) + 2 ] • ( 4 - 6 ) 

The subtraction term, on the other hand, diverges logarithmically for m —? 0. 
As shown in appendix D the divergence is such that the limit 

P2 = { ( 4^)2 + ( 4J ) 

exists. Numerically one obtains 

P 2 = 0.02401318111946489(1). (4.8) 

Putting these results together, we conclude that eq.(4.2) holds except perhaps 
for the form of the error term (at this point we only know it vanishes for 
P o). 

Eq.(4.2) may also easily be derived using coordinate space methods. The 
starting point is the representation 

A{p) = lim Y e~ ex2 e~ ipx G{x ) 2 . (4.9) 

xeA 

We then note that the function H(x) [eq.(2.7)] satisfies 

1(4 x 4 

H(x) ~- 2 ( 2^2 \ 1 - 2 + 8 T~tA 

x^oo TT Z (X Z ) Z [ X z (x z ) J 

7 x 4 x 6 (x 4 ) 2 ) 

“ Jtff +i0 ~(Yy ~ 288+28 °-^Jy + •••) ( 4 ‘ 10 ) 

and so has the same leading order behaviour as G(x) 2 . The convergence of 
the sum (4.9) can thus be accelerated by separating a singular part A s (p) 
according to 

A(p) = A s (p) +A r (p), (4.11) 
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where 


A s (p) = -(4tt)- 2 lim V e~ ex2 e~ ipx H (x), (4.12) 

e—^ 

A r (p) = y e-^ [G(x f + (4ir)- 2 H(x)\ . (4.13) 

xeA 

Since A s (p) is just the Fourier transform of H(x) and since the square bracket 
in the last expression above falls off like | a; | —6 , we deduce that 

A (p) p z 0 -yy ln (p 2 ) + A Ao) + 0(p 2 ). (4.i4) 

We have thus rederived eq.(4.2) with the constant P 2 being given by 

^2 = ~7^V2 + 22 [G(x) 2 + (dvr) 2 H(x)] . (4.15) 

Moreover the argument shows that the error term is as quoted. The next-to- 
leading terms in the expansion could in fact easily be calculated by making 
a further subtraction. We finally remark that the sum in eq.(4.15) can be 
computed using the techniques of sect. 3. 


4-2 The bigmac diagram 

We now proceed to discuss the integral 

d 4 k d 4 g 


D{p) = J 


1 


/_ T (2vr) 4 (2vr) 4 k 2 g 2 r 2 ' 

(cf. fig. lb). Our goal is to establish the expansion 


r = —k — g — p 


(4.16) 


D(p) = A 0 +p 2 [A t + Bt ln(p 2 )] + 0(p 4 ) 

p —^0 


(4.17) 


and to compute the coefficients A\ and B\ (the leading term Aq has already 
been calculated in sect. 3). 

The situation is very much the same as in the case of the one-loop integral 
A(p) discussed above. We first note that 

D(p) = y e- vpx G{x f (4.18) 

xeA 
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and split the sum into a singular and a regular part according to 

D{p) = D s {p) + Drip), (4-19) 

where 

Dsip) = -|(4vr)- 4 Y, e”^A H(x), (4.20) 

xeA 

D r {p) = Yj i G ^ 3 + H 4lr )~ 4AH ( X )\ ■ (4.21) 

xeA 

The splitting has been done in such a way that the square bracket in the last 
expression falls off like |s| _8 at large x [cf. eqs.(3.12),(4.10)]. The expansion 
of D r {p) 1° order p 2 is thus obtained straightforwardly by expanding the ex¬ 
ponential factor in eq.(4.21). As for the singular part we note that D s ip) is 
just the Fourier transform of AH(x), i.e. D s ip) is proportional to p 2 ln(p 2 ). 
Taken together our discussion shows that the expansion (4.17) holds with 

A 1 = x2 [ G ( x f + H 4ir )~ 4AH ( x )\ (4-22) 

xeA 

and B\ = 4(47t) -4 . To calculate the sum (4.22) we again make use of the 
techniques introduced in sect. 3. Note that Hix) is computable at all non¬ 
zero x through eq.(2.6). 77(0) is calculated in appendix D. In this way the 
result 

Ai = -0.00007447695(1) (4.23) 

is obtained. 


4-3 The eye diagram 

The integral corresponding to the scalar eye diagram (fig. lc) may be written 
in the form 

, , r d 4 k i 

E(p) = Lw¥?- Mk) ’ (4 ' 24) 

where Aik) denotes the one-loop integral studied previously and s = k + p as 
before. In the following our aim is to show that 


Dip) = 


o 2(4vr) 


[!n(p 2 )]‘ 


1 

(47t) 2 



ln(p 2 ) + X t + 0(p 2 ) (4.25) 
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and to calculate the constant X\. 

We begin by defining the function 


A as (k) 


1 

(47r) 2 


ln(£; 2 ) + 2 


P 2 - 


As discussed in subsect. 4.1 the one-loop integral A(k) satisfies 


A(k)-A as (k) k = Q 0(k 2 ). 


(4.26) 


(4.27) 


This suggests that we introduce singular and regular parts of the eye diagram 
through 


^ , . 1 r d 4 k In (k 2 ) 

sip) -~X^f k 2 P ’ 


(4.28) 

(4.29) 


The relation 


E(p) = E s (p ) + E r (p) + 



A(p) 


(4.30) 


is an immediate consequence of these definitions. 

We first consider the regular part E r (p). The function A(k) — A as (k ) is 
smooth in the whole integration range except at k = 0 where its behaviour is 
given by eq.(4.27). Taking this into account it is easy to show that 


E r (p) = E r (0) + 0(p 2 ). 

p —*■(_) 


(4.31) 


The constant E r ( 0) can be computed using position space techniques. From 
the discussion in subsect. 4.1 we infer that 


A(k) - A as (k) = Y, ( e ~ tkX ~ !) [G( x ? + (4 tt)- 2 H(x)\ . (4.32) 

xeA 

As already noted before the square bracket in this expression falls off like 
\x\~ & a t large x. When eq.(4.32) is inserted in eq.(4.29), the integration and 
summation may hence be interchanged and it follows that 

E r { 0) = K(x) [G{x) 2 + (4t t)~ 2 H(x)\ . (4.33) 

xeA 
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The auxiliary function K(x) has been introduced in subsect. 3.3 and is dis¬ 
cussed further in appendix C. At large x the terms in the sum (4.33) have an 
asymptotic expansion consisting of two series of the usual form, of which one 
is multiplied by ln(a; 2 ) [cf. eq.(C.8)]. The subtraction method of sect. 3 may 
be applied in such cases as well. We only have to introduce the derivatives of 
the appropriate zeta functions to calculate the subtraction term. The result 

E r { 0) = -0.000002546129(1) (4.34) 


is thus obtained. 

We still need to evaluate the singular part E s (p). This is a one-loop inte¬ 
gral which may be calculated using the momentum space subtraction technique 
previously employed in the case of the integral A(p) (subsect. 4.1). We thus 
split the integral into three parts, 


Es(p) 


{E' S (P, m) + E”(m) + In (m 2 )A(p)} , 


where 


E' s {Pi m) 


j"* d 4 k j ln(k 2 /m 2 ) ln(k 2 /m 2 + 1) | 

(2vr) 4 1 p]2 (jfc2 +m 2)2 J ’ 


(4.35) 


(4.36) 


E"( 


rri 


j"* d 4 k In (k 2 /m 2 + 1) 
-x (2vr) 4 (jfc2 + m 2)2 


(4.37) 


The auxiliary mass parameter m > 0 is again scaled to zero in a fixed propor¬ 
tion to p. From the Reisz power counting theorem we then infer that 


lim 

A^O 


E' s (Xp,Xm) 


00 d 4 k fin (k 2 /m 2 ) ln(k 2 /m 2 + 1) 1 
(27t) 4 \ k 2 s 2 ( k 2 + m 2 ) 2 J 


= - { \ [M™ 2 Ip 2 )] 2 + ln(m 2 /p 2 )} . (4.38) 


The second part, E"(m ), is logarithmically divergent for m 0 in such a way 
that the limit 


Ps = 


lim { - - —- 

■o 1 2(4vr) 


[ln(m 2 ) 


P 2 


(4vr) 


In (m 2 ) + E”(m)> (4.39) 
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exists. This is shown in appendix D where P 3 is also computed numerically. 
Collecting all contributions to eq.(4.35) the result 


E °W A t ln(p2) ] 2 “ ln(p2) } “ i4A P3 + 0(p2) 


( 4*0 


(4.40) 


is thus obtained. 

We hnally go back to eq.(4.30) and conclude that the expansion (4.25) 
holds with 


X! 


E r (0) + 




(4.41) 


Numerically one obtains 


Xi = 0.001123187953(1). (4.42) 


As a last point we mention that X\ has previously been computed by Alles 
et al. [5] through numerical integration in momentum space. Their result was 
estimated to be accurate to 4 significant decimal places but differs from the 
number given above by 10%. We do not know the cause for this discrepancy. 


5. Concluding remarks 


In this paper attention has been restricted to lattice theories in four di¬ 
mensions with free propagators given by eq.(1.3). The extension of our meth¬ 
ods to other situations is however straightforward. In particular, the recursive 
computation of free propagators, such as the Dirac propagator in lattice QCD 
with Wilson quarks, is generally possible along the lines discussed in sect. 2. 
An amusing remark in this connection is that in two dimensions the propagator 
of a scalar massless held with standard lattice action is at any point x exactly 
computable as a linear combination of 1 / 7 r and 1 with rational coefficients. 

Most diagrams that we encountered in our computation of the relation 
between the bare lattice coupling and the MS coupling in SU(iV) gauge theo¬ 
ries [1] have integrands with complicated numerators. It turns out, however, 
that the techniques introduced in this paper are sufficient to compute all of 
them. The important point to note is that the continuum limit of most of 
these integrals is rather easy to take. The resulting expansion coefficients are 
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integrals of the type discussed in sect. 3 and so can be calculated by the meth¬ 
ods described there. Further details on these computations will be published 
elsewhere. 

We are indebted to Claus Vohwinkel for sharing his insights about the 
free propagator with us. 


Appendix A 


In this appendix we derive the large x expansion of the free propaga¬ 
tor given in eq.(2.4). The methods employed are more generally useful. We 
begin with the expansion of some auxiliary functions and then consider the 
propagator. 

A.l Auxiliary functions 

Let Q(p) be a homogeneous polynomial in the momentum p with degree d, 
and n > 1 an integer such that d — 2n > — 3. The integral 

/ oo 14 

Pe^e- p 2 Q(p)( P 2 )- n (A.l) 

-OO W 

is then absolutely convergent and well-defined for all x E 1R 4 . Our aim in the 
following is to establish the behaviour of I(x) for large x. 

Consider first the case n = 1. Here we have 

r °° r 00 r| 4 w . 2 

'<*)=«(-»)/, d< L(w e ' p ' e " ,p 

= Q(-id)(4ir 2 x 2 )- 1 ^1 - e~ x , (A.2) 

from which it follows that 

I(x) ~ <Q(-i<9)(47r 2 a; 2 ) -1 (A.3) 

x —OO 

up to exponentially small corrections. 
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Let us now assume that n > 2. Define to = 2n — 4 and let T m be the 
Taylor operator of order to, viz. 

m ^ 

T m f(x) = A 0 )- ( A - 4 ) 

k =0 

Since d > to we have 

Q(p)e ipx = Q(-id)( 1 - T m )e ipx (A.5) 


and so conclude that 

1 Z’ 00 

I(x) = Q(-id)-, -T-T / d t(t- 1)"- 1 (4tt/)- 2 (1 - T m )e- 2 / 4t . (A.6) 

Note that the differentiation and integration could be interchanged here be¬ 
cause (1 — T m )e tpx e~ p ( p 2 )~ n is absolutely integrable. Next we substitute 
/—►!// and use the binomial formula. This yields 


I(x) 


1 \n-l~n 


X 


f dtt~ p |e- ta;2 / 4 - ^^(-^ 2 /4)4 • (A.7) 

4° l i/=o z/ ' J 


Through partial integration all terms with fj, > 1 can be reduced to 

f dtt~ k (e~ tx / 4 — l) ~ — ln(a; 2 ) + constant + ... 

Jo \ / ^-xx> 

In this way we obtain the asymptotic expansion 
I(x) 


1 f t-l)™- 1 4 


(n — 1)! x 2 


o ( 47 t ) 2 

n — 1 

5 (P~ l) ! d!( ra “ 1 -d)! V 4 

/4— -L 


(- 1 ) 


71 — 1 


2 \ Ai — 1 


ln(a; 2 ) >, (A.8) 


which is valid up to exponentially small corrections. Note that after differ¬ 
entiation, the logarithm disappears and one is left with terms of the form 
P(x)(x 2 )~ k , where P is a polynomial and k an integer. 
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A.2 Large x expansion of the propagator 
Let h(p), p E 1R 4 , be a smooth function with 


h(p) 


1 if |p| < 1, 
0 if \p\ > 2. 


(A.9) 


G{x) then coincides with 


/ oo i 4 


to all orders of I#I . Next we note that 


(A.fO) 


f = E( _1 ) 


71+1 


71—1 

where we have introduced the notation 

3 


(2 ra )! 


P 2n , 


(A.lf) 


p n = £ov) n - 

/ i =0 


(A-12) 


It follows that 


e p \p 2 ) 1 p = Q {p 2 ) 1 + 1 + P 4 (p 2 ) 2 + tp 2 + p 4 (p 2 ) 1 

- p 6 (p 2 )~ 2 H-(p 4 ) 2 (p 2 ) -3 + • • • 

360 ^ J 144 J J 


(A.13) 


When this series is multiplied with e p and subtracted from the integrand in 
eq.(A.lO), the singularity at p = 0 is smoothed. As a result we have 


G(x) 


d 4 p 


1 1 p 4 


~ *°° ./-oo ( 27r ) 4 Ip 2 12 (p 2 ) 2 

and hence, using the results of subsection A.l, 

G(l) ™ ++ - +'‘ lll(I 4 


(A.14) 


-^d & \n{x 2 )P ^^{d^fx 2 ^ 2 )) P0{\x\ 8 ). (A.15) 
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Carrying out the differentiations one obtains the expansion (2.4). Higher order 
terms could be generated easily, using a symbolic manipulation program. 


Appendix B 


In this appendix we establish the relations (2.12) and (2.13) following an 
argument first given by Vohwinkel [2]. The starting point is the observation 
that the recursion relation (2.9) becomes one-dimensional along the lattice 
axes. Explicitly, if we define 


9o(n) 

= G(n, 0, 0,0), 

(B.l) 

9i(n) 

= G(n, 1, 0,0), 

(B.2) 

02 «) 

= G(n, 1,1,0), 

(B.3) 

0 a(ra) 

= G(n, 1,1,1), 

(B.4) 

and use the lattice symmetries, it 

is easy to show that 


9o(n + 1) = 8g 0 (n ) - 60 !(n) - g 0 (n - 1), 

(B.5) 

2?? 

9i( n + 1) = n+1 [ 40 i(n) ■ 

l 

tQ 

0 

(5) 

1 

to 

tQ 

to 

(5) 

1 

s s 
+ 1 

1 

(B. 6 ) 


2t % ji — 2 

02 (n + 1) = [402 (rc) - 201 (n) - g 3 (n)\ - ^-p-^0 2 (ra - 1), (B.7) 

O77 77 Q 

g 3 (n + 1) = —— [403(n) - 30 2 (n)]-— ~g 3 {n - 1), (B.8) 

n + 3 n + 3 


for all n > 1. These equations look deceptively simple, but so far we have not 
been able to solve them exactly. Vohwinkel has however noticed that 

h(n) = (n - l)go(n) + Sng^n) + 3 (n + l) 0 2 (n) + (n + 2 )g 3 (n) 

- ng 0 (n - 1) - 3(n - l) 0 r(n - 1) 

- 3(n - 2 )g 2 (n - 1) - (n - 3 )g 3 (n - 1) (B.9) 
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is a “constant of motion”, i.e. using the recursion relations one may show that 


ki(n + 1) = ki(n) for n > 1. (B.10) 

Another such constant of motion is given by 

k 2 (n) = n 2 ki(n) — 2 (n — 1 )(n — 2)go(n) — 3n(2n — 3)</i(n) 

- 6 (n + l)(n - l)g 2 {n) - (n + 2)(2 n - 1 )g?,{n) - 3 ng 0 (n - 1) 


- 6(n - l)gi(n - 1) - 3 (n - 2)g 2 (n - 1), (B.ll) 

but no further independent constant was found. 

The values of k\ and k 2 are easily worked out in the limit n oo where 

g^in) = (47r 2 ra 2 ) -1 + 0(l/ra 3 ) for all g = 0,..., 3. (B.12) 

The result is 

h = 0, (B.13) 

k 2 = -4/7r 2 , (B.14) 


and setting n = 1 in eqs.(B.9) and (B.ll) now directly yields eqs.(2.12) and 
(2.13). 


Appendix C 


C.l Properties of K(x) 

From the definition (3.26) of the auxiliary function K(x) it is straightforward 
to show that 

-A K{x) = G{x) (C.l) 

and 

(Vl + V^)K(x) = -x^G(x). (C.2) 
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Using the latter equation we can determine K(x) at any point x recursively 
in terms of G and the values of K at the corners of the unit hypercube at the 
origin. These are given by 


Tt (0, 0,0,0) = 0, (C.3) 

id(l,0,0,0)= — Ig(0,0,0, 0), (C.4) 

id(l, 1,0,0)= —^G(0, 0,0,0), (C.5) 

id(l,1,1,0)= -^G(0, 0,0,0)+ 0,0), (C.6) 

A'(l, 1,1,1)= ^(0,0,0,0)- 0,0)- (C.7) 


K(x) is thus obtained as a linear combination of the form (2.14) and so is 
accurately computable. 

At large x the asymptotic expansion 

K{x) ™ { hl| - T2) + c = + 0(| - T| "‘ 1 } |C ' 81 

holds, where 

ci = 3.560146537656400(2). (C.9) 

Eq.(C.8) has been derived by applying the techniques discussed in appendix 
A. The leading singularity at p = 0 in the defining integral (3.26) (which gives 
rise to the logarithm in the large x expansion) was taken care of by introducing 
an appropriate auxiliary function I(x). 


C.2 Properties of L(x) 

The function L(x) defined through eq.(3.29) satisfies 

-AL(x) = K(x) + ±G(0)-^ 


and 


(V* + V^)L(x) 


1 

— x u 
2 ^ 


K(x) + ±G(0) 


(C.10) 


(C.ll) 
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Its values at the corners of the unit hypercube at the origin are 


£(0,0,0,0)= 0, (C.12) 

£(1,0,0,0) = -^G(0,0,0,0)+ 2 ^ 2 , (C.13) 

£(1,1,0,0)= --^G(0, 0,0,0)+-^, (C.14) 
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£ ( 1 , 1 , 1 ,0 ) =- —G(0,0,0, 0 ) +^2, (C.15) 


£(1,1,1,1)= -^G(0,0,0,0)+ ^G(l, 1,0,°)+ 1^2- ( C ' 16 ) 

This function is, therefore, also expressible as a linear combination of the form 
(2.14). 

The asymptotic behaviour of L(x) is given by 

L ( x ) ,=00 ^ - ln (® 2 ) + c 3 + + 0(|m|- 2 )| 



(C.17) 

c 2 = -0.49811600254393(2), 

(C.18) 

c 3 = -2.77499211742106(1). 

(C.19) 


Appendix D 

The computation of the constants P 2 , P 3 and H( 0) discussed in this 
appendix is based on the identity 

/ 7r 1 4 u -1 r co 

v (2^ (P + m2) '“ = T Xajl [e- 2t /o(2f)] 4 , (D.l) 


24 



which holds for all a > 0 and m > 0. The function Iq{z) appearing in this 
formula is a modified Bessel function (ref.[6], §8.43). To evaluate the constant 


P 2 = lim 

771^0 


(4tt) 


■ ln(m 2 


d 4 k 1 

(2vr) 4 (jfc2 +m 2)2 


(D.2) 


we insert eq.(D.l) with a = 2 and discuss the limiting behaviour of the integral 
for m 0. The Bessel function Iq{z) is a smooth function with asymptotic 
behaviour such that 


e“ 2t / 0 (2t)l 4 = (47rt)“ 2 + 0(t“ 3 ). (D.3) 

L J t —* OO 


It follows from this that the integral diverges as m 0, but it is now easy to 
show that the limit (D.2) exists and is given by 

p 2 = + J o dtt { [e- 2t I 0 (2t)} 4 - 0(t - I)(47r7)- 2 } , (D.4) 


where 7 = —r'(l) denotes Euler’s constant. From this representation it is 
straightforward to compute P 2 by numerical integration, using library routines 
for the Bessel function. As a result one obtains the number quoted in eq.(4.8). 

In the case of the constant P3 [eq.(4.39)] we differentiate the identity 
(D.l) with respect to a at a = 2 and proceed as above. In a few steps, taking 
eq.(D.4) into account, one obtains 


Ps = 


(47r) 


:7 “7 


7r" 

12 


f°° 

+ J dtt[l - 7 - ln(/)] j [e _2t / 0 (2f)] 4 - 0{t - l)(47i7 ) -2 j , (D.5) 
and numerical integration now yields 


P 3 = 0.03467243216643385(1) (D.6) 


(for better convergence the series expansions of the Bessel function have been 
used in the range 0 < t < 0.01 and 100 < t < 00). 

We finally discuss the computation of the constant 


H( 0) = 


d 4 k 


In (k 2 


(D.7) 
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In this case we first restrict the parameter a in eq.(D.l) to the range 0 < a < 1. 
We can then set m = 0 in the formula. After performing a subtraction at / = 0 
the identity becomes 



1 

r(a + 1) 


1 f°° 

+ f{ajl dUa_1 { [ e “ 2t/o(2 ^ 4 “ 9(1 ~ *)} ’ (D ' 8) 

which is now valid for — 1 < a < 1. Differenting at a = 0 we obtain 

f°° 

#(0) = -7-y d //" 1 {[e- 2t / 0 (2f)] 4 - d(l-t)} (D.9) 

and numerical integration gives 


H{ 0) = 1.99970764451731256(1). (D.10) 
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